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ABSTRACT 

We consider the statistical relationship between the growth rate of activity 
in the early phase of a solar cycle with its subsequent amplitude on the basis 
of four datasets of global activity indices (Wolf sunspot number, group sunspot 
number, sunspot area, and 10.7-cm radio flux). In all cases, a significant corre- 
lation is found: stronger cycles tend to rise faster. Owing to the overlapping of 
sunspot cycles, this correlation leads to an amplitude-dependent shift of the solar 
minimum epoch. We show that this effect explains the correlations underlying 
various so-called precursor methods for the prediction of solar cycle amplitudes 
and also affects the prediction tool of Dikpati et al. (2006) based upon a dynamo 
model. Inferences as to the nature of the solar dynamo mechanism resulting from 
predictive schemes which (directly or indirectly) use the timing of solar minima 
should therefore be treated with caution. 

Subject headings: SumActivity 



1. Introduction 

Solar activity is the driver of space weather, which has practical consequences for human 
activities in space. This is one reason why the search for methods to predict its (short-term 
and long-term) future levels has found much interest in the literature. Further motivation 
arises from the potential implications for understanding the origin of solar activity: a reliable 
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method of predicting the amplitude of future solar cycles could provide a constraint on 
dynamo models. On the other hand, the con verse proposition is not necessarily true, as has 
been pointed out by lBushby fc Tobias! (120071 ): the nonlinear dynamics of the dynamo might 
be such as to make mid- to long-term prediction impossible even if an almost perfect physical 
understanding of the dynamo mechanism is achieved. 

In many cases, recipes for prediction are inferred from correlations found in historical 
records of measu red quantities, which are (directly or i ndirectly) related t o sola r activity. 
As illustrated by Iwilsonl (11994 cf. his Fig. 14.2) and lLantos & Richard! (119981 , cf. their 
Fig. 6), the success of most methods in actually predicting the unknown amplitude of a future 
cycle is rather disappointing. Nevertheless, the correlations between several 'precursors', i.e., 
quantities measured during the descending or minimum phase of a cycle and the amplitude of 



the subsequent cycle (e.g., lHathaway et al.lll999l ; ISchattenll2003l ) might have a non-random 



origin and thus call for a physical explanation. This could have implications for dynamo 
models. 

In this paper, we consider the effect of the overlapping of solar cycles in combination 
with their asymmetric shape on the correlations between precursors and following cycle 
amplitudes. In this connection, the important aspect of the asymmetry is the difference 
of the amplitude-dependent ascent rate in the early cycle phase (related to the so-called 
Waldmeier effect) compared to the decay rate n ear the end of a cycle. Since sunspot cycles 
overlap for typically 2 to 3 years (iHarveyi Il992l ). this asymmetry affects the timing of the 
activity minima, which are pivotal epochs for most precursor methods. We show that these 
effects can explain the correlations upon which such methods are based, without necessarily 
implying a direct physical connection between the precursor quantity and the following cycle. 
We also show that the essence of the Waldmeier effect, i.e., that stronger cycles tend to show 
a faster rise of activity levels during their ascending phase than weaker cycles, is a robust 
property present in all activity indices. On this basis, we explain how cy cle asymmetry and 
cycle o verlapping may also affect the dynamo-ba sed prediction method of iDikpati &: Gilman 
( 120061 ). in spite of recent claims to the contrary (IDikpati et al.ll2008bl ). 



2. Overlapping asymmetric cycles and precursors 

Precursor methods are based upon the existence of a correlation between some physical 
quantity measured during the descending or minimum phase of a cycle and the amplitude of 
the following cycle. The correlation is established by considering historical records of data, 
the longest of which is the record of sunspot numbers. Assuming a non-random origin of 
such a correlation, there are three possible explanations for its existence: 
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The precursor is a feature of the old cycle that represents or is related to an input 
quantity for the dynamo process, so that its magnitude directly affects the amplitude 
of the next cycle. For example, it has been suggested that the strength of the polar 
field during solar minimum is a measure of the poloidal field from which the toroida l 
field for the next cycle is generated (e.g.. ISchatten et al.lll978l ; IChoudhuri et al.ll2007l ). 



The precursor quantity represents an early manifestation of the new cycle, which al- 
ready affects the high latitudes of the S un but does not ye t produce sunspots. The 
notion of such an 'extend ed solar cycle' (IWilson et al.lll988l ) is support ed by observa - 
tions of coronal activity (lAltrock Il997l). ephemeral m agnetic regions (jHarveyl Il994l ). 



and zonal flows ( Howe et al. 



2006 



Altrock et all 120081 ). 



If the definition of the precursor quantity directly or indirectly depends on the timing 
of the solar minimum, the amplitude-dependent shift of the minimum due to the over- 
lap of asymmetri c cycles can lead to a correla tion of the precursor with the amplitude 
of the new cycle ( ICameron fc Schusslerl 120071 . see their Fig. 10). This does not neces- 
sarily involve a physical connection between the precursor and the dynamo process, in 
principle even permitting the prediction of random cycle amplitudes. 



The third possibility arises from the observation that informati on about the height of the 
maximum of a cycle is already contained in its early rise phase (jWaldmeierl Il935l . Il955l ). In 
fact, the shapes of most historical solar cycle s can reasonably well b e d escrib ed by simple 



functions containing a few parameters (e.g., lHathaway et al.l Il994 ; iLil Il999l ). so that an 



estimate of the fu r ther development of a cyc le is possible a few years after sunspot minimum 
( jWaldmeierl 1 19361 ; lElling fc Schwentekl Il992l ). Since cycles overlap, i.e., the rise of the new 
cycle starts already when the decay of the old cycle is still ongoing, the timing of the cycle 
minimum (the epoch when the sum of the activities of both cycles is minimal) and its 
height are affected by the amplitude of the new cycle: the minimum occurs earlier and is 
higher for a subsequent strong cycle than for a weak cycle. In fac t, stronger cycles tend 
to be preceded by shorter cycles (earlier minima, ISolanki et al.ll2002l ) with higher minimum 
activity levels. Both quantities can be considered as precursors with corr elation coefficients 
with the amplitude of the next cycle of about 0.7 (lHathaway et al.lll999n . 



Cameron fc Schusslerl (120071 ) have noted that the amplitude-dependent shift of the mini- 



mum epoch is enhanced by the amplitude- dependent asymmetry of solar cycles: strong cycles 
tend to grow faster in activity than weak cycles, while the decay rate in the late descending 
phase is largely independent of the cycle amplitu de. This is related to, but not identical to, 
the so-called 'Waldmeier effect' ( I Waldmeierl 1 1 9 351 ) which is often stated in th e form: the time 



between cycle minimum and maximum is shorter for stronger cycles (e.g., lHathaway et al. 
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20021 ). However, the crucial quantity for the minimum shift is not the time between mini- 



mum and maximum but the steepness of the activity rise (decay) in the initial (late) phase 
of the cycle, which can be determined independently of the timing of the activity minima 
and maxima by considering the slope of the activity curve. We show in the next section 
that the correlation between initial growth rate and cycle amplitude is a very robust feature 
shown by all available global activity indices. 



3. The robust amplitude-dependent cycle asymmetry 



Hathaway et al. found that the 'classical' Waldmeier effect, i.e., the correlation 



between the rise time (time interval between minimum and maximum) an d the cycle ampli 
tude, is by a factor of about two weaker for the group sunspot numbers (IHoyt fc Schatten 
19981 ) than for the Wolf sunspot numbers. For both datasets, the historically strongest cycles 
(numbers 18, 19, 21, and 22) do not fit the linear relationship very well, showing somewhat 
long rise times. This weakening of the correlation can be understood by the shift in the 
timing of the minima: stronger, rapidly growing cycles lead to early preceding minima, thus 
increasing the time between minimum and ma ximum. This effect be comes even more pro- 
nounced in the analysis of sunspot area data by lDikpati et al.l (l2008bl ) owing to the fact that 
the sun spot areas for cyc les 21, 22 and 23 reach a slightly higher second 'Gnevyshev max- 
imum' (iGnevyshevi 119671 ) a few years after the first maximum coinciding with the sunspot 
number maximum. It is not surprising that determining the rise time in terms of these later 
maxima destroys the correlation with the cycle amplitude. 

As we have pointed out in the previous section, the relevant quantity for the minimum 
shift of overlapping cycles is not the time interval between minimum and maximum but 
the steepness of the activity growth ( decay ) in the initial ( late) phase of the cy cle. An 
inspection of Fig. 5 of lHathaway et al.l (120021 ) and of Fig. 1 of iDikpati et al.l (l2008bl ) already 
suggests that the different datasets do not differ strongly in this respect. In order to give 
a quantitative account, we need to measure the steepness without reference to the epochs 
of maximum and minimum; in particu lar, the l atter is already affected by the correlations 
that we wish to study. For instance, iLantod (120001 ) considered the slope of the sunspot 
number curve at the inflexion point during the ascending part of cycles 9 to 22 and found 
a correlation coefficient of r = 0.88 with the cycle amplitude, while the (anti) correlation 
with the rise time from minimum to maximum only yields \r\ = 0.61. In our case, we need 
to consider specifically the early rise and late decay phases, so that we estimate the rise 
and decay rates by determining the time required for the activity index under consideration 
to cover a fixed interval of values. This interval is chosen such that it is not significantly 
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affected by the cycle overlap. The rise and decay rates are then defined as the ratios of the 
value intervals and the corresponding rise and decay times, respectively. 

We have analyzed four datasets: monthly Wolf sunspot numberj j] (SN) for cycles 7 
to 23, monthly group sunspot number s ^ (GS N, Hoyt fc Schatten 19981 ) for cycles 7 to 22, 



sunspot areas (SAR, iBalmaceda et al.l 120051 ) for cycles 12 to 23, and 10.7-cm solar radio 



flu^ (SRF) for cycles 19 to 23. While we have used the full available data sets for SRF 
and SAR, the SN and GSN data are only considered from cycle 7 on because a) we wish 
to compare the results for the different datasets during roughly the same period of time, 
and b) for very low-amplitude cycles like those during the Dalton minimum and the first 
cycles after the Maunder minimum we cannot use the same intervals for the determination 
of the rise and decay times. Fig. [T] shows the datasets, also indicating the intervals chosen for 
the determination of growth and decay rates: 30-50 (SN, GSN), 300-600 microhemispheres 
(SAR), and 900-1100 sfu (SRF). All data have been smoothed by a Gaussian with a FWHM 
of 1 year. 



1 http:/ /side. oma.be 

2 http:/ /www. ngdc.noaa.gov/stp/SOLAR/ftpsunspotnumber. htm 
3 http: / /www.drao.nrc.ca/icarus/www/sol_home.shtml 
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Fig. 1. — Data sets used in this study: Wolf sunspot number (SN, top left), group sunspot 
number (GSN, top right), sunspot area (SAR, bottom left), and 10.7-cm solar radio flux 
(SRF, bottom right). SN and GSN are considered from cycle no. 7 onward (after the Dalton 
minimum). The dashed horizontal lines indicate the intervals used for the determination of 
the rise and decay rates. 
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Fig. 2. — Scatter diagrams of rise rates (left panels) and decay rates (right panels) versus 
cycle amplitude for four sets of solar activity indices: Wolf sunspot number (SN, top row), 
group sunspot number (GSN, second row), sunspot area (SAR, third row, cycles 21-23 
indicated by triangles, two of which nearly coincide on the left-hand panel), and 10.7-cm 
solar radio flux (bottom row). The corresponding linear correlation coefficients are given. 
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Scatter diagrams of the resulting rise and decay rates with the cycle amplitudes are 
shown in Fig. [21 together with the corresponding correlation coefficients. The figure demon- 
strates that high correlation coefficients (r = 0.83 . . . 0.89) between rise rate and cycle am- 
plitude exist for all activity indices. Assuming a Gaussian distribution, the corresponding 
significance levels are between 95% (SRF) and over 99% (SN, GSN, SAR). In particular, the 
group sunspot numbers and the sunspot areas exhibit an even higher correlation than the 
Wolf sunspot numbers. Note that cycles 21-23 (indicated by triangles in the third row of 
Fig. [2]), which, owing to their higher Gnevyshev maxima in the SAR data, largely destroye d 



the correlation of the 'classical' Waldmeier effect in the analysis of iDikpati et al.l (l2008bl ). 
nicely fit the correlation for the growth rate. On the other hand, the decay rates show 
no significant correlation with the cycle amplitude, so that we find an amplitude- dependent 
asymmetry between the rise rates early in a cycle and the decay rates in the late phase. 
This asymmetry acts to make the behavior near sunspot minimum rather insensitive to the 
amplitude of the old cycle but sensitive to that of the new. 

We have checked how strongly the correlation coefficients depend on the chosen intervals 
for the determination of the growth and decay rates, and found them to be rather insensitive. 
For instance, in the case of the growth (decay) rates for the SN data, we have r = 0.79 (0.05) 
for the interval [25,45], r = 0.82(0.22) for [30,50] (the interval chosen for Fig. [2]), r = 
0.85 (0.42) for [35, 55], and r = 0.88 (0.54) for [40, 60]. Similar results are found for the other 
datasets. 



4. Implications for precursor methods 

The overlapping of cycles with an amplitude-correlated asymmetry leads to a depen- 
dence of the minimum epoch on the amplitude of the following cycle: since the initial activity 
of a large-amplitude cycle increases faster while the decline in the decay phase is not or less 
amplitude-dependent, the minimum level of activity between overlapping cycles occurs the 
earlier the higher the amplitude of the following cycle. Such a minimum shift affects the 
value of any precursor quantity which is defined (directly or indirectly) with regard to the 
minimu m epoch; in the extreme, it m ay completely explain the predictive power of the pre- 



cursor. 



Cameron fc Schiisslerl (120071 ) have demonstrated this for the case of a very simply 
defined precursor, namely, the activity level (such as sunspot number or area), three years 
before minimum. The earlier the minimum occurs (i.e., the stronger the new cycles is), the 
higher is the precursor value, simply because it is measured at an earlier time when the 
activity of the old cycle had not declined as much. This explains why the method can, in 
principle, even predict cycles with random amplitudes. It is important to note, however, that 
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such methods always require that the minimum epoch is already known. For a real prediction 
(as opposed to a postdiction of past cycles) such a method can be applied only some time 
into the new cycle, when proper averaging can be performed to define the minimum tim e in 
view of the strong fluctuations of activity around minimum (cf. lHarvey fc White! Il999l ). 



Precursors that are defined (directly or indirectly) in relation to the timing of the min- 
imum are affected by the amplitude-dependent minimum shift, which thus may (partly or 
fully) explain their correlation with the amplitude of the next cycle (i.e., their predictive 
power). Among the various precursor methods affected we mention just a few examples: geo 



magnetic activity before and around solar minimum (e.g.. lOhllll966l ; lLegrand &: Simonlll 981 
Lantos fc Richard! Il998l ; iKand 120071 ) , sunspot activity ar ound minimum (IHathaway et al. 
999 ), length of th e preceding su nspot cycle (jKand 120081 ) . skewness of the preceding cycle 



cycle (IPodladchikova et al 



dRamaswamyj|l977l; iLantoa 120061) . as well as properties of the decay phase of the preceding 



2008 



5. Implications for the model of Dikpati et al. 



Cameron fc Schiisslerl (120071) s uggested that at least pa rt of the predictive skill of the 
model of iDikpati et al.l (120061 ) and iDikpati fc Gilmanl (120061 ) . henceforth referred to as the 
DGT model, may result from the overlap of asymmetric cycles in the sunspot area data from 
which their source term is derived. We have shown above that the underlying effect, i.e., 
the amplitude-dependent asymmetry between the growth and decay phases, is present in all 
global activity datasets, including the sunspot area data used in the DGT model. 

The effect of this crosstalk between cycles on the DGT model results from the assump- 
tion of a fixed speed and range of latitude drift of the source term in the course of a cycle. 
The shortening of the cycle length (defined as the time between minima) by the shift of 
the minimum preceding a strong cycle then leads to larger source amplitudes (sunspot ar- 
eas) assigned to the low-latitude part of the artificial butterfly diagram used in the DGT 
model. As a consequence, the magnetic flux crossing the equator as the potentially relevant 
quantity for the dynamo amplitude of the next cycle is higher, while the opposite is true 
for a weak following cycle. In this way, the correlation in the data survives the artificial 
stretching/ compressi ng of cycles carried ou t in the DGT model in order to obtain a constant 
cycle period. In fact. IDikpati et al.l (j2008al ) found that, in their model, the flux crossing the 
equator during a cycle is significantly correl ated (r = 0.76) with the amplitude of the next 



cycle. This is in accordance with Fig. 8a in IDikpati &: Gilmanl (120061 ). which demonstrates 
that the memory in th eir model does not ex tend much longer than one cycle, but it is in 
contrast to the claim of IDikpati et al.l (j2008al ) that it takes 17-21 years for the poloidal sur- 
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face field to reach the tachocline and generate the new toroidal field, so that 2-3 preceding 
cycles contribute to the toroidal flux of a cycle. A stronger diffusive coupling between 
surface and tachocline than indicated by the explicit value of their magnetic 
diffusivity (for e xample, due to n umerical diffusion) could possibly explain this 
discrepancy (cf. lYeates et ajj 12008). 



By disregarding the observed sunspot latitudes in favor of an imposed fixed latitude 
drift, the DGT model becomes affected by the amplitude-dependent minimum shift. The 
extent to which this contributes to the predictive skill of the model could be easily tested 
by replacing the actual sunspot area data as input by a synthetic dataset with random cycle 
amplitudes (without any m emory), but keeping t h e corr elation resulting from the minimum 
shift of overlapping cycles. ICameron k, Schiisslerl (120071 ) used this approach to demonstrate 
the effect on precursor methods. The amount of predictive skill shown by the DGT model 
for such an input will directly indicate how strongly this correlation affects the prediction. 



6. Discussion 

The shift of the solar minimum epoch due to the amplitude-dependent asymmetry be- 
tween the growth and decay phases of solar cycles can explain the rather high correlation 
coeffi cients between various precursor quantities and the strength of the subsequent cycle 



^e.g., lHathaway et al.lll999l ). In this connection, we should note that even a statistically 
highly significant correlation does not automatically imply a high skill of the precursor for 
the actual prediction of a future cycle: 



1. A typical value of, say, r = 0.8 for a 'good' precursor leaves about one third of the 
variance in amplitude unexplained by the correlation, so that the prediction is prone to 
co nsiderable s t atisti cal uncertainty. Even the value of r = 0.97 shown by the method 
of iThompsonl (119931 . see Hathaway et al. 1999) did not prevent an utterly inaccurate 
prediction for cycle 23 (30% too high). 

2. The data used to determine the correlations in most cases cover only part of the known 
range of variability in the solar cycle. Periods such as the Maunder and Dalton minima 
indicate that nonlinear and stochastic effects may severely limit the predictability since 
it is unclear how st rong ly these effects are represented in the observed surface flux 



(IBushby fc Tobias! 120071 ). 'Tuning' of free parameters in prediction schemes can lead 
to seemingly high correlations for past cycles while the predictive skill for future cycles 
in fact is much lower. 
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The dependence of many precursor methods on knowing the minimum epoch further 
limits their practical usefulness: the large relative fluctuations of the activity indices around 
sunspot minimum make it necessary to average the data (over 1 year, for instance) in order 
to obtain a proper definition of the minimum epoch. Consequently, a sensible prediction 
can only be made one or two years after the minimum has passed, so that there is no big 
advantage compared to methods that use the informati on provided by the early phases of the 
new cycle itself, like the steepness of the activity rise (j Waldmeierl 1 1 9 36l ; lEUing fc Schwentek 
1992h . 



The relation between the shift of the minimum time and the amplitude of the next 
cycle due to cycle overlapping is strongest for a dynamo where the time interval between the 
onsets of activity of subsequent cycles is constant. Fluctuations in the timing of the onset 
of activity introduces noise into the relation between minimum shift and activity and thus 
reduces the predictive skill of precursor methods. Quantifying such fluctuations in the solar 
cycle data is difficult because the length of th e available da t asets is not s ufficient to reliably 
determine the phase fluctuations and drifts ( Goughl 1978 ; Hovng Il996l ) and to decide on 



the existence or otherwise of a solar 'clock' (cf. lDickdll978l . Il988l ). However the expected 
shifts due to the overlapping of the cycles and the amplitude dependent asymmetry produce 
variations in the cycle length (as measured from minimum to minimum) which are consistent 
with those observed since the Dalton minimum. This means the effect we have outlined is 
relevant. In passing, we also note that in the case of advection-dominated dynamo models 
t he cycle period is determined by the meridional flow, w hich supports phase stability 
jCharbonneau &: Dikpatil I2OO0I : Icharbonneaul [2003 ) . 



For our understanding of the solar dynamo, it would certainly be useful to clarify to 
which extent the skill of some prediction methods results 1) from correlations in the input 
data (such as discussed here), 2) from capturing early high-latitude manifestations of an 
extended new cycle, and 3) from capturing mechanisms that connect properties of the old 
cycle to the strength of the new cycle (such as a dynamo model with a memory of at least one 
cycle). It is reasonable to suppose that there is some memory in the solar cycle i.e., that the 
solar cycle amplitudes do not constitute a purely random sequence. Potentially important 
information could be gleaned from identifying a quantity that unequivocally represents this 
memory and is not 'contaminated' by early information leaking in from the new cycle. 
Clarifying whether the memory extends over more than one cycle w ould help to decid e 
between diffusion-dominated and advection-dominated dynamo models (jYeates et al.ll2008l ). 
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7. Conclusions 

We have confirmed a highly significant correlation between the growth rate of activity 
during the early phase of a solar cycle and its maximum amplitude for all global activity 
indices, i.e., Wolf and group sunspot numbers, total sunspot area, and 10.7-cm radio flux. 
On the other hand, there is no significant correlation between the decay rate in the late cycle 
phase and the cycle amplitude. 

Owing to the overlapping of individual cycles, this asymmetry leads to an amplitude- 
dependent shift of the minimum epoch, thus explaining (fully or partly) the predictive power 
of precursor methods which (directly or indirectly) use the timing of the activity minimum 
as a pivotal point. The resulting correlation in t he sunspot area da t a pro bably also affects 



the predictions with the dynamo-based model of iDikpati fc Gilmanl (120061 ) 



For our understanding of the origin of the solar magnetic field, it is important to disen- 
tangle the effects of 'real' physical precursors, i.e., properties of the old cycle directly affect- 
ing the flux generation for the next cycle or early high-latitude manifestations of the new 
cycle, from apparent precursors, which derive their predictive power from the the amplitude- 
dependent shift of the minimum epoch. 
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